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FOREWORD 


A  mission  of  the  National  Imagery  and  Mapping  Agency  (NIMA)  requires  that  geodetic 
quality  absolute  positions  be  determined  autonomously  and/or  in  areas  far  removed  from  other  sites. 
The  use  of  Differential  Global  Positioning  systems  (DGPS)  techniques  or  relative  positioning  is 
sometimes  precluded.  Keyed  GPS  receivers  can  output  Precise  Positioning  Service  (PPS)  navigation 
solutions.  Because  they  are  computed  in  real  time,  PPS  solutions  necessarily  use  the  broadcast 
ephemerides  and,  consequently,  are  subject  to  ephemeris  and  satellite-clock  prediction  errors.  The 
PPS  solutions  can  be  recorded,  but  the  broadcast  ephemerides  and  the  observations,  when  corrected 
for  Selective  Availability  (SA),  cannot.  This  report  demonstrates,  through  simulations,  that  field- 
recorded  PPS  solutions  can  be  reprocessed  with  the  more  accurate  precise  ephemerides.  Reprocessing 
with  the  precise  ephemerides  offers  improved  navigation  solutions,  but  does  not  require  that  the  SA- 
corrected  pseudorange  or  phase  observations  be  saved.  The  method  is  called  Precise  Absolute 
Navigation  (PAN).  This  report  continues  the  evaluation  of  the  PAN  algorithm,  which  has  been 
described  in  previous  documentation. 

Support  for  this  development  was  provided  by  NIMA  under  the  direction  of  Mr.  Stephen 
Malys  (NIMA/TRB)  and  Mr.  Ken  Croisant  (NIMA/TRB,  ESMD). 

The  author  wishes  to  thank  Dr.  Randall  Smith  (NIMA/TR)  for  suggesting  this  technique 
as  one  of  interest  to  NIMA  for  particular  applications.  The  author  also  wishes  to  thank  the 
following  individuals  at  NIMA  in  St.  Louis  for  assisting  in  the  development  of  PAN:  Mr.  Gus 
Jones,  Mr.  Brian  Tallman  (GIMGH),  and  Mr.  Brian  Deem,  for  providing  data  and  suggestions; 
and  Ms.  Lisa  McCormick  for  processing  data  sets  through  the  Data  Correction  Facility  (DCF). 

This  report  has  been  reviewed  by  Mr.  Everett  Swift,  Program  Manager  for  NIMA-related 
projects;  Dr.  Jeffery  Blanton,  Head,  Space  Systems  Applications  Branch;  and  Mr.  James  Sloop, 
Head,  Space  and  Weapons  Systems  Analysis  Division. 


Approved  by; 


IARD  T.  LEE,  Acting  Head 
Theater  Warfare  Systems  Department 
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1.0  INTRODUCTION 

The  purpose  of  the  Precise  Absolute  Navigation  (PAN)  development  is  to  improve,  by 
postprocessing,  the  accuracy  of  a  set  of  Precise  Positioning  Service  (PPS)  navigation  solutions  from 
keyed  receivers.  Real-time,  stand-alone  navigation  solutions  necessarily  use  the  broadcast 
ephemerides  as  the  source  of  the  satellite  position  and  satellite  clock  estimates.  However,  in  certain 
cases,  the  particular  mission  could  benefit  from  more  accurate  navigation  solutions.  If  the  mission  is 
not  time  critical,  it  may  be  acceptable  to  wait  until  the  postfit  ephemerides  are  available  and,  thus, 
improve  the  PPS  navigation  solutions  with  these  precise  ephemerides.  In  this  case,  significantly  better 
results  can  be  expected. 

With  precise  ephemerides,  static  absolute  position  solutions  can  achieve  submeter  repeatability 
[1],  [2].  Meter-level  navigation  solutions  have  been  demonstrated  by  postprocessing  observations 
with  postfit  precise  ephemerides,  clock  estimates,  and  an  ionospheric  model  [3].  Applications  where 
precise  absolute  navigations  would  be  of  interest  include  cases  where  auxiliary  data  for  differential 
or  relative  solutions  are  not  available,  or  where  the  track  of  a  moving  vehicle  or  ship  is  desired  at 
remote  locations  for  from  other  sites. 

In  the  operational  environment,  the  PPS  navigation  solutions  will  be  obtained  from  a  keyed 
receiver.  A  keyed  receiver  will  be  able  to  correct  for  the  effects  of  Selective  Availability  (SA)  and 
AntiSpoofing  (AS)  that  are  intentionally  introduced  to  degrade  the  accuracy  of  the  system.  AS 
encrypts,  and  SA  corrupts  the  Global  Positioning  System  (GPS)  signals— -making  the  precise  signal 
and  the  full  accuracy  of  the  broadcast  ephemerides  accessible  to  authorized  users  only.  These  keyed 
receivers  can  display  the  navigation  solution  in  real  time,  but  the  field-corrected  satellite  ephemerides 
and  the  SA-corrected  pseudorange  and  phase  observations  are  classified  and  therefore  are  not 
available  to  the  user.  The  PPS  position  solutions  are  unclassified  and  can  be  recorded  in  the  field, 
along  with  the  GPS  time  and  the  particular  satellites  contributing  to  the  solutioa  In  order  to  improve 
the  PPS  solutions  at  a  later  time  with  the  postfit  precise  ephemerides,  the  following  information  is 
required: 

>•  the  GPS  epoch 

>-  the  satellites  contributing  to  the  PPS  solution  at  that  epoch 

>  the  PPS  solution 

>•  the  broadcast  ephemerides  used  to  form  the  PPS  solutions 

>-  the  precise  ephemerides 

Section  2  begins  by  comparing  PAN  solutions  with  three  relative  positioning  methods  using 
three  data  sets  that  explore  various  problems  that  could  be  experienced  with  real-world  conditions. 
Wherever  possible,  On-The-Fly  kinematic  relative  positioning  is  used  as  truth  for  these  evaluations. 
The  accuracy  of  kinematic  positioning  has  been  shown  to  be  in  the  centimeter  level.  This  has  been 
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documented  elsewhere,  as  seen  in  References  [4],  [5],  and  [6].  Section  3  continues  with  an  analysis 
of  a  four-day  data  set  collected  by  the  National  Imagery  and  Mapping  Agency  (NIMA)  over  a  road 
course  between  St.  Louis,  Missouri,  and  Toledo,  Ohio,  in  July  1997.  No  reference  sites  dedicated  to 
this  experiment  were  set  out,  but  some  data  from  the  National  Oceanographic  and  Atmospheric 
Administration  (NOAA)  Continuously  Operating  Reference  Sites  (CORS)  and  the  Naval  Surface 
Warfare  Center,  Dahlgren  Division  (NSWCDD),  Dahlgren,  Virginia,  were  used  to  obtain  relative 
positioning  solutions. 


2.0  PAN  RESULTS  COMPARED  WITH  OTHER  METHODS 

This  section  investigates  how  well  the  absolute  terrain  heights  along  selected  ground  tracks  can 
be  determined.  Previous  work  has  shown  that  in  the  case  of  airborne  vehicles,  PAN  is  capable  of 
±l-m  determination  of  ellipsoid  heights.  PPS  navigation  solutions,  as  well  as  relative  methods,  are 
also  possible  candidates  for  the  proposed  task. 

Depending  upon  the  accuracy  required  for  the  terrain  heights,  any  of  five  different  positioning 
methods  could  be  used.  PPS  navigation  solutions  and  PAN  navigation  solutions  derived  from  the  PPS 
are  the  two  methods  that  do  not  require  a  reference  site.  Three  relative  positioning  methods  that 
require  at  least  one  static  reference  site  are  range-corrected  Standard  Position  Service  (SPS),  relative 
positioning  with  one  or  two  frequency  pseudoranges,  and  kinematic  relative  positioning.  The 
advantage  of  the  relative  methods  is  that  all  data  can  be  handled  unclassified  and  collected  without 
a  keyed  receiver. 

To  perform  an  evaluation,  two  existing  data  sets  were  identified  as  test  cases,  and  a  third  data 
set  was  collected  locally  near  Naval  Surface  Warfare  Center,  Dahlpen,  Virginia  (NSWCDL).  The 
two  existing  data  sets  include  aircraft  data  collected  by  the  Applied  Research  Laboratory  at  the 
University  of  Texas  at  Austin  (ARL/UT)  in  1995,  which  has  been  used  in  previous  evaluations  of 
PAN,  and  the  data  set  that  was  collected  at  Holloman  Air  Force  Base  (AFB)  by  NSWCDD,  NIMA 
and  several  other  organizations  inNovember  1 996.  The  static  sites  that  collected  data  simultaneously 
with  Holloman  include  the  Naval  Postgraduate  School  (NPS)  Monterey,  ARL/UT  at  Austin, 
NSWCDD  at  Dahlgren,  and  four  locations  in  Canada. 

The  third  data  set  was  collected  by  placing  an  antenna  on  a  vehicle  roof  and  traveling  along 
Route  301  near  Dahlgren.  Two  40-km,  round-trip  circuits  of  this  course  were  made  to  evaluate 
repeatability  of  the  solutions. 


2.1  Piper  Aircraft  Data  from  ARL/UT 

The  Piper  aircraft  data  was  collected  by  ARL/UT  on  day  237  of  1 995.  It  is  a  well-behaved  data 
set  with  observations  collected  at  5-s  intervals.  For  this  analysis,  603  epochs  were  used  to  perform 
the  comparison.  The  On-The-Fly  kinematic  relative  positioning  technique  was  used  as  the  truth.  The 
time  span  of  the  analysis  was  limited  by  the  kinematic  integer  ambiguity  resolution  restrictions  caused 
by  the  distance  between  the  aircraft  and  the  reference  site,  as  illustrated  in  Figure  1 .  The  aircraft  flew 
over  the  reference  site  MSI  en  route  to  its  loiter  area.  This  allowed  the  integers  to  be  established. 
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Figure  1 — Piper  Aircraft  Ground  Track 


However,  the  aircraft  distance  was 
mostly  greater  than  30  km  from  MSI 
and,  as  a  consequence,  if  satellite  signals 
dropped  out  or  were  lost,  the  integers 
could  not  be  reestablished.  The  east- 
north  ground  track  of  the  Piper  aircraft 
is  shown  in  Figure  1  with  respect  to 
MSI,  and  the  ellipsoid  height  is  shown 
in  Figure  2.  The  height  plot  indicates 
that  the  altitude  was  continually  chang¬ 
ing  on  the  order  of  10  to  20  m,  probably 
due  to  turbulence.  This  effect  can  simu¬ 
late  the  height  variations  of  terrain  and 
so  is  a  good  test  case  for  that  reason. 
The  geometric  dilution  of  precision 
(GDOP)  and  the  number  of  satellites  in 
view  at  the  aircraft  are  plotted  in 
Figure  3. 


Position  differences  were  com¬ 
puted  among  the  four  alternative  pro¬ 
cessing  methods  and  the  On-The-Fly 
kinematic  solutions,  which  were  taken 
as  truth.  The  four  methods  were: 


( 1 )  Relative  positioning  using  smoothed 
pseudoranges 

(2)  Simulated  PPS  navigation  solutions 

(3)  PAN-improved  PPS  solutions  using 
the  precise  ephemerides 

(4)  Range-corrected  SPS  navigation 
solutions. 


Cases  (1),  (4),  and  (5)  (Truth)  require 
the  use  of  a  reference  site.  The  other 
two  do  not  need  a  reference.  However,  PPS  needs  SA-corrected  data  (a  keyed  receiver);  PAN  needs 
the  PPS  solutions,  a  list  of  the  satellites  used  at  each  epoch,  a  copy  of  the  broadcast  ephemerides  used 
by  the  receiver  to  compute  the  PPS  solutions,  and  a  precise  ephemeris. 

The  differences  in  the  east,  north,  and  vertical  directions  are  plotted  in  Figures  4  (PPS),  5 
(PAN),  6  (range-corrected  SPS),  and  7  (smoothed  pseudorange  relative  positioning).  Notice  that 
Figures  4  and  6  have  a  maximum  vertical  scale  of ±15  m,  while  the  other  two  have  a  scale  of  ±5  m. 
Based  on  the  standard  deviations,  the  smallest  deviations  from  the  truth  come  from  the  smoothed 
pseudorange  solutions,  followed  by  PAN,  range-corrected  SPS,  and  PPS.  Dr.  Benjamin  Remondi’s 


PIPER  AIRCRAFT  DATA,  DAY  237,  1995 
EAST-HEIGHT  TRUTH  GROUND  TRACK 


Figure  2 — Piper  Aircraft  Ellipsoid  Height 
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smoothed  relative  pseudorange  solu¬ 
tions  are  a  good  choice  if  keyed  receiv¬ 
ers  are  to  be  avoided.  In  this  case,  the 
reference  site  can  be  thousands  of  kilo¬ 
meters  distant.  This  long-distance  capa¬ 
bility  will  be  demonstrated  in 
Section  2.4. 

This  data  set  illustrates  that  under 
normal  conditions  where  satellite  track¬ 
ing  is  steady  and  GDOP  is  good,  the 
solutions  from  all  techniques  are  accept¬ 
able.  In  particular,  the  PAN  algorithm 
did  a  good  job  of  improving  the  PPS 
solutions.  Comparison  of  the  standard 
deviations  listed  in  Figures  4  and  5 
shows  a  reduction  in  each  component 
by  almost  a  factor  of  four. 


- - - N 

PIPER  AIRCRAFT  DATA,  DAY  237,  1995 
SATELLITE  TRACKING  AT  AIRCRAFT 


l  _ J 

Figure  3 — Piper  Aircraft  GDOP  and  the  Number  of 
Satellites  Being  Tracked 
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PIPER  AIRCRAFT  DATA,  DAY  237,  1995 
SIMULATED  PPS  MINUS  TRUTH 


Figure  4 — PPS  Solutions  Minus  Truth  for  the  Piper  Aircraft 
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PIPER  AIRCRAFT  DATA,  DAY  237,  1995 
PAN  IMPROVEMENT  MINUS  TRUTH 


Figure  5— PAN  Solutions  Minus  Truth  for  the  Piper  Aircraft 


PIPER  AIRCRAFT  DATA,  DAY  237,  1995 
RANG  E  CO  RRECTED  SPS  MINUS  TRUTH 


Figure  6 — Range-Corrected  SPS  Minus  Kinematic  Truth  for  the 


Piper  Aircraft 
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PIPER  AIRCRAFT  DATA,  DAY  237,  1995 
SMOOTHED  PSEUDORANG  ES  MINUS  TRUTH 


GPS  TIME  OF  WEEK  (s) 

_ J 

Figure  7 — Smoothed  Pseudorange  Relative  Positioning  Minus  Kinematic 

Truth  for  the  Piper  Aircraft 


2.2  Dahlgren/Port  Royal  Route  301 

Circuits 

A  map  of  the  Dahlgren  region  is  shown  in 
Figure  8.  A  vehicle  with  an  Ashtech  Z-XII 
receiver  and  a  GPS  antenna  on  its  roof  collected 
1-s  data  on  1  April  1997.  A  reference  site  was 
operating  simultaneously  at  NSWCDL.  The 
route  the  vehicle  took  begins  at  NSWCDL  and 
continues  to  Route  301  and  northeast  to  the 
Potomac  River.  At  the  Potomac,  the  vehicle 
reversed  course  and  traveled  south  along  Route 
301  to  Port  Royal.  At  the  crossroads,  the  vehicle 
turned  around  and  headed  back  north  to  the 
Potomac.  This  route  was  repeated  a  second 
time. 

In  contrast  to  the  aircraft  data,  the  Route 
301  data  was  far  from  benign;  it  is  now  question¬ 
able  whether  it  was  suitable  for  evaluating  the 
merits  of  the  various  processing  schemes.  There 
were  many  satellite  dropouts  due  to  blockages 


Figure  8 — Map  of  Dahlgren  and  Port  Royal 
Area 
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by  passing  vehicles  and  roadside  trees  and  terrain.  The  blockages  resulted  in  repeated  losses  of  lock 
of  satellites  by  the  Ashtech  Z-XII  receiver.  This  caused  havoc  with  the  processing  programs,  which 
previously  had  not  been  stressed  to  this  degree.  The  lesson  learned  is  that  data  collected  along  a 
highway  route  may  be  extremely  difficult  to  process  unless  a  concerted  effort  is  expended  to  build 
adequate  protection  into  the  existing  software.  Getting  meaningful  operational  results  from  this  type 
of  data  may  require  considerable  hand  editing,  given  the  current  maturity  of  the  software. 

The  GDOP  and  the  number 
of  satellites  tracked  during  the 
two  circuits  of  the  course  are 
plotted  in  Figure  9.  Compare  this 
figure  with  its  counterpart  for  the 
aircraft  solutions  in  Figure  3.  The 
repeated  loss  and  reacquisition  of 
satellites,  as  shown  by  the  dashed 
line  with  x  symbols  in  the  figure, 
were  mainly  due  to  obstructions 
by  trees  along  the  roadside.  Sev¬ 
eral  of  the  satellites  available  dur¬ 
ing  this  time  were  at  low  eleva¬ 
tion  angles,  and  at  certain  places 
the  trees  significantly  obscured 
the  sky.  The  GDOP,  indicated  by 
the  solid  line,  is  dependent  upon 
the  number  of  satellites  available 
at  any  time  and  fluctuates  as  the 
number  in  view  changes.  This 
causes  the  accuracy  of  the  solutions  to  change  from  epoch  to  epoch  as  the  satellites  come  in  and  out 
of  view.  During  the  first  half  of  the  data  span,  the  number  of  satellites  in  view  rarely  surpassed  five, 
but  there  were  fewer  dropouts  than  later  in  the  span.  During  the  second  half  of  the  span,  the  number 
of  satellites  visible  in  the  sky  was  greater,  but  the  number  of  dropouts  was  more  frequent.  Eventually, 
the  On-The-Fly  kinematic  truth  could  not  obtain  correct  integers  for  the  new  satellites  and  stopped 
producing  solutions.  Figure  10  shows  the  gaps  in  the  truth  during  the  second  circuit.  This  figure 
shows  the  changing  distance  between  the  vehicle  antenna  and  the  static  reference  site  located  on  the 
roof  of  Building  1200  at  NSWCDL. 

Since  the  interest  here  is  in  the  accuracy  of  the  determination  of  ellipsoid  height,  Figure  1 1 
shows  the  height  solutions  obtained  from  each  of  the  five  positioning  techniques  plotted  on  the  same 
graph.  The  plot  covers  the  entire  time  period  of  observations.  The  kinematic  solutions,  which  are 
considered  the  most  accurate  when  they  exist,  are  shown  as  the  filled  circles.  At  the  beginning  of  the 
observations  (at  about  229,500  s),  the  integer  ambiguities  have  not  been  established,  and  the 
kinematic  ellipsoid  heights  are  at  zero  (no  solutions).  After  a  few  hundred  seconds,  the  kinematic 
solutions  begin  and  continue  uninterrupted  until  about  235,200  s.  There  are  then  two  periods  of 
dropouts  as  indicated  in  Figure  10.  The  other  four  positioning  techniques  produce  solutions  as  long 


GDOP  AND  NUMBER  OF  SATELLITES  TRACKED 
PPS  SOLUTIONS,  RT-301:  1  APRIL  1997 


Figure  9 — GDOP  and  the  Number  of  Satellites  Being  Tracked 

by  the  Receiver 
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as  there  are  enough  unblocked 
satellites  in  view.  In  a  very  gen¬ 
eral  sense,  all  the  solutions  agree 
to  within  about  10  m.  When  there 
are  bad  solutions,  such  as  near 
233,600  s,  they  are  obviously  bad. 

Repeatability  can  be  dem¬ 
onstrated  by  viewing  a  short  seg¬ 
ment  of  the  route  on  both  of  the 
two  circuits.  On  Route  301,  the 
vehicle  covered  the  same  ground 
in  the  same  lane  each  time 
around,  so  the  heights  are  repro¬ 
ducible.  The  only  differences  were 
in  the  geometry  of  the  satellites, 
the  speed  of  the  vehicle,  and  stops 
due  to  traffic  signals.  Figure  12 
shows  a  segment  of  the  south¬ 
bound  course,  which  can  be  iden¬ 
tified  in  Figure  1 1  by  comparing  the  time  of  week.  Figure  13  shows  the  same  segment  on  the  second 
lap.  Particular  features  due  to  variations  in  the  road  height  can  be  identified  on  both  plots. 
Qualitatively,  the  repeatability  for  the  kinematic  solutions  looks  good  to  within  1  -  m  peak  difference. 
This  gives  some  credibility  for  using  the  kinematic  as  truth  under  these  conditions.  Presumably,  the 
error  is  as  large  as  it  is  because  of  the  constant  changing  of  the  satellites  in  view.  It  would  be 
interesting  to  investigate  in  detail  the  behavior  of  the  pseudorange  and  phase  observations  at  the 
vehicle,  but  this  would  take  a  considerable  amount  of  time  to  do.  With  a  better  editing  algorithm 
available,  it  may  be  possible  to  eliminate  many  of  the  spikes  seen  in  the  solutions.  They  are  probably 
due  to  the  existence  of  questionable  observations  making  their  way  into  the  solutions. 

Differences  among  the  four  techniques  and  the  kinematic  truth  are  shown  in  Figures  14  through 
1 7.  All  are  plotted  as  east,  north,  and  vertical  components  on  a  ±1 0-m  scale.  The  solution  spikes  are 
due  to  bad  observations  as  described  above.  The  simulated  PPS  and  the  PAN  solutions  are  similar 
in  character,  with  the  PAN  solutions  being  slightly  less  noisy  as  determined  by  the  standard 
deviations.  In  this  data  set,  the  ephemerides  are  not  the  major  source  of  error.  The  PAN  algorithm 
uses  the  PPS  solutions  as  the  starting  point  and  exchanges  the  broadcast  ephemerides  for  the  precise. 
Any  errors  in  the  PPS  solutions  not  due  to  the  ephemerides  are  carried  unchanged  over  into  the  PAN 
solutions.  The  gap  in  all  plots  between  235,200  and  236,000  s  is  where  the  truth  is  lacking. 


-  . 

KINEMATIC  RANGE  FROM  REFERENCE  SITE 
RT-301  EXCURSION:  1  APRIL  1997 
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Figure  10 — The  Range  of  the  Vehicle  from  the  Reference  Site 
Showing  the  Periods  of  Dropouts 


2.3  Positioning  the  NIMA  Vehicle  at  Holloman  AFB 

During  the  evening  of  7  November  and  the  early  morning  hours  of  8  November  1996,  two 
NIMA  vehicles  carried  Ashtech  receivers  to  various  known  sites  at  Holloman  AFB  in  New  Mexico. 
A  NIMA  reference  site  was  also  operating,  but,  due  to  a  software  problem,  the  data  from  this  receiver 
was  not  saved.  However,  during  about  2  of  the  12  hours  of  NIMA  operation,  NSWCDD  operated 
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ELLIPSOID  HEIGHTS 


RT-301  EXCURSION:  1  APRIL  1997 


Figure  11 — Ellipsoid  Height  Solutions 


ELLIPSOID  HEIGHTS,  FIRST  CIRCUIT 
RT-301  EXCURSION:  1  APRIL  1997 


Figure  12 — Ellipsoid  Heights  for  a  Southbound  Segment  on  Lap  1 
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ELLIPSOID  HEIGHTS,  SECOND  CIRCUIT 
RT-301  EXCURSION:  1  APRIL  1997 


Figure  13 — Ellipsoid  Heights  for  a  Southbound  Segment  on  Lap  2 


SIMULATED  PPS  SOLUTIONS  MINUS  KINEMATIC  SOLUTIONS 
RT-301  EXCURSION:  1  APRIL  1997 


Figure  14 — PPS  Solutions  Minus  Kinematic  Truth 
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PAN  SOLUTIONS  MINUS  KINEMATIC  SOLUTIONS 
RT-301  EXCURSION:  1  APRIL  1997 


GPS  TIME  OF  WEEK  (s) 

Figure  15 — PAN  Solutions  Minus  Kinematic  Truth 


ERECTED  SPS  SOLUTIONS  MINUS  KINEMATIC  SOLUTIONS 
RT-301  EXCURSION:  1  APRIL  1997 


Figure  16 — Range-Corrected  SPS  Minus  Kinematic  Truth 


NSWCDD/TR-98/61 


r - - - a 

DIFFERENTIAL  SOLUTIONS  MINUS  KINEMATIC  SOLUTIONS 
RT-301  EXCURSION:  1  APRIL  1997 


V _ _ _ / 

Figure  17 — Smoothed  Pseudorange  Relative  Positioning  Minus  Kinematic 

Truth 


a  receiver  at  site  F901  not  far  from  the  NIMA  vehicle  route.  This  2-hr  segment  was  used  here  to 
evaluate  the  ellipsoid  height  accuracy  at  a  closed  course  and  at  previously  surveyed  markers. 


A  plot  showing  the  route  of 
the  NIMA  vehicle  that  was  on  the 
course  during  the  2-hr  period  is 
shown  in  Figure  18.  The  route  be¬ 
gins  at  the  upper  left,  and  progresses 
east  and  south  toward  the  lower 
right,  stopping  at  each  of  the  sites 
indicated.  Near  site  ST11,  the 
vehicle  turned  around  and  headed 
north  again.  The  2-hr  period  ended 
just  north  of  ST63.  The  NSWCDD 
reference  site  F901  is  shown  as  a  ® 
symbol  at  the  origin  near  the  lower 
right  comer  of  the  plot.  The  GDOP 
and  the  number  of  satellites  in  view 
is  shown  in  Figure  1 9.  At  this  time  of 
day  and  at  this  location,  the  number 
of  satellites  tracked  at  F901  were 
never  fewer  than  seven,  and  the 


Figure  18 — NIMA  Vehicle  Route 
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GDOP  was  three  or  smaller  most  of 
the  time. 


This  data  set  provides  an  op¬ 
portunity  to  compare  results  against 
independently  established  absolute 
positions.  The  markers  indicated  in 
Figure  18  have  known  coordinates. 
Two  of  these  will  be  used  to  illus¬ 
trate  the  accuracy  and  repeatability 
of  the  kinematic  relative  positioning 
technique.  Of  course,  any  error  in 
the  position  of  the  reference  site 
F901  will  add  to  the  apparent  abso¬ 
lute  error  in  the  other  site  positions. 


The  NIMA  operation  involved 
carrying  a  GPS  antenna  attached  to 
the  outside  of  the  vehicle,  so  that 
data  was  collected  continuously  while  in  motion  as  well  as  when  stopped  at  a  marker.  The  antenna 
on  its  pole  would  be  detached  from  the  vehicle  and  walked  to  the  marker.  The  point  at  the  bottom 
of  the  pole  would  be  placed  on  the  monument  center  and  the  pole  leveled.  It  would  stay  stationary 
there  for  several  minutes.  It  would  then  be  walked  back  to  the  vehicle  and  reattached.  The  vehicle 
would  then  proceed  on  to  the  next  site,  where  the  process  would  be  repeated. 

An  example  of  the  process 
described  above  is  shown  in  Fig¬ 
ure  20  at  site  PBM2.  This  plot 
shows  the  east  and  north  compo¬ 
nents  of  the  stop  at  PBM2  with 
respect  to  F901  in  meters.  The 
known  location  of  the  mark  is 
indicated  as  a  *  symbol  in  the 
plot.  The  kinematic  relative  posi¬ 
tion  of  the  marker  is  at  the  cluster 
of  dots  about  0.6-m  west  and 
0.3-m  south  of  the  cross.  Abso¬ 
lute  errors  in  the  tabulated  loca¬ 
tion  of  PBM2  and  F901  probably 
account  for  the  disagreement  in 
position.  Figure  21  shows  the  east 
and  ellipsoid  height  components. 
Again,  the  PBM2  mark  is  indi¬ 
cated  by  the  *  symbol  and  the 
antenna  location  (2.067  m  above 
the  mark)  is  indicated  by  the  tri- 
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WADGPSAT  HOLLOMAN  AFB,  11/8/96 
A  PORTION  OF  THE  NIMA  TRUCK  ROUTE  -  SITE  PBM2  - 
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Figure  20 — NIMA  East-North  Antenna  Position  in  the 
Vicinity  of  Site  PBM2 


Figure  19 — GDOP  and  the  Number  of  Satellites  in  View  at 

Site  F901 


13 


ELLIPSOID  HEIGHT  (m) 


NSWCDD/TR-98/61 


angle.  The  agreement  in  height  is 
very  good.  Figure  22  shows  a 
high-resolution  look  at  the  scatter 
in  the  kinematic  solution  in  the 
east  and  ellipsoid  height  compo¬ 
nents.  A  1-cm  diameter  circle 
(plotted  as  an  ellipse  because  of 
the  different  scales  on  the  two 
axes)  is  included  for  reference. 
The  repeatability  of  the  kinematic 
(truth)  solutions  during  the  time 
the  antenna  was  over  the  mark  is 
at  the  centimeter  level  as  claimed. 

The  stops  at  site  ST63  will 
be  displayed  next.  This  is  the  only 
site  during  the  2  horns  of  F901 
operation  where  two  stops  were 
made.  Figure  23  shows  the  east- 
north  activity  in  the  vicinity  of 
ST63,  and  Figure  24  shows  the 


Figure  22 — Scatter  in  the  East-Fleight  Position  of  PBM2  from  Kinematic 

Positioning 
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north-ellipsoid  height  components.  The  first  stop  was  in  the  southbound  direction  and  the  second  was 
in  the  northbound  direction.  The  accepted  position  for  ST63  is  indicated  by  the  ©  symbol.  The  actual 
position  of  the  mark  appears  to  be  near  the  horizontal  blob  in  the  center  of  Figure  23. 

A  closer  look  at  the  solutions  while  the  antenna  was  at  the  mark  is  shown  in  Figures  25  and  26. 
The  southbound  solutions  are  shown  as  filled  circles  and  the  northbound  solutions  as  open  circles. 
A  3-cm  diameter  circle  is  included  on  each  plot  for  reference.  The  filled  symbols  form  a  tighter  cluster 
within  the  reference  circle  than  the  open  symbols,  but  both  sets  are  included  within  the  circle.  This 
shows  that  the  kinematic  truth  solutions  are  repeatable  at  the  3-cm  level  after  a  period  of  time  has 
passed,  and  the  constellation  of  satellites  have  formed  a  different  configuration. 

The  evidence  just  presented  for  sites  PBM2  and  ST63  shows  that  the  kinematic  truth  is  reliable, 
so  it  will  be  used  to  evaluate  the  positioning  accuracy  of  the  other  four  techniques  during  the  2-hr 
period.  Figures  27  through  30  show  the  differences  between  the  simulated  PPS  solutions,  PAN,  range 
corrected  SPS,  and  smoothed  pseudorange  solutions,  respectively;  and  the  truth.  These  differences 
are  smaller  than  those  from  the  aircraft  flight  and  all  are  plotted  in  ±5-m  axes. 

The  solutions  minus  the  truth  from  the  three  test  data  sets  discussed  in  Section  2  are 
summarized  in  Tables  1  through  4.  The  standard  deviation  in  each  component  from  the  three  test 
cases  is  computed  and  listed  in  the  Average  a  row  of  each  table.  The  Route  301  standard  deviations 
are  much  greater  than  the  others  and  dominate,  but  are  consistent  in  relative  magnitude  with,  the 
other  two.  If  the  Route  301  results  are  discarded,  all  four  techniques  have  mean  values  smaller  than 
10  m  from  the  truth  and  standard  deviations  smaller  than  5  m  The  PAN  means  are  1  m  or  smaller, 


Figure  23 — Kinematic  Solutions  in  the  Vicinity  of  ST63  in  the  East-North 

Components 
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Figure  24 — Kinematic  Solutions  in  the  Vicinity  of  ST63  in  the  East-Height 

Components 


WADGPS  AT  HOLLOMAN  AFB,  11/8/96 
A  PORTION  OF  THE  NIMA  TRUCK  ROUTE  -  SITE  63 
, _ KINEMATIC  POSTIONING  WITH  RESPECT  TO  SITE  F901  . 


SYMBOLS  EACH  SECOND 


— —  SOUTHBOUND 
— •—  NORTHBOUND 


3  cm  RADIUS  CIRCLE 


NORTH  (m) 


Figure  25 — Scatter  in  the  East-North  Solutions  at  Site  ST63 
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WADG  PS  AT  HOLLO  MAN  AFB,  11/8/96 
A  PORTION  OF  THE  NIMA  TRUCK  ROUTE  -  SITE  63 


NORTH  (m) 


Figure  26 — Scatter  in  the  North-Height  Solutions  at  Site  ST63 


SIMULATED  PPS  MINUS  KINEMATIC  TRUTH 
ST6WITH  RESPECT  TO  F901 


GPS  TIME  OF  WEEK  (s) 


Figure  27 — PPS  Solutions  for  the  NIMA  Vehicle  Minus  Truth 
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Figure  28 — PAN  Solutions  for  the  N1MA  Vehicle  Minus  Truth 


RANGE  CORRECTED  SPS  MINUS  KINEMATIC  TRUTH 
ST6WITH  RESPECT  TO  F901 


GPS  TIME  OF  WEEK  (s) 


Figure  29 — Range-Corrected  SPS  Solutions  Minus  Truth 
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Table  2 — PAN  Improvement 


Units  = 

EAST  1 

NORTH 

VER1 

_l 

< 

O 

meters 

MAX 

MIN 

MEAN 

a 

MAX 

MIN 

MEAN 

a 

MAX 

MEAN 

o 

PIPER 

0.90 

-0.75 

0.05 

0.42 

1.82 

-0.40 

0.54 

0.52 

2.89 

-0.60 

1.02 

0.83 

RTE  301 

-257.99 

108.15 

-2.06 

6.70 

264.64 

-104.35 

1.42 

6.87 

158.51 

-116.13 

-3.19 

5.06 

HOLLOMAN 

0.61 

-2.01 

-0.82 

0.35 

0.22 

-1.78 

-0.59 

0.40 

0.95 

-4.35 

-0.76 

0.81 

AVERAGE  a 

2.49 

2.60 

i 

2.23 
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Table  3 —  Range-Corrected  SPS 


Units  = 

EAST 

NORTH 

VERTICAL 

meters 

MAX 

MIN 

MEAN 

a 

MAX 

MEAN 

a 

MAX 

MIN 

MEAN 

a 

PIPER 

2.17 

-0.58 

0.48 

0.52 

0.64 

-3.32 

-1.19 

0.65 

5.81 

-5.64 

0.90 

1.95 

RTE  301 

243.60 

-255.40 

-0.49 

14.74 

488.2 

-249.3 

1.46 

24.17 

189.7 

-118.6 

-0.03 

17.57 

HOLLOMAN 

1.88 

-14.00 

0.03 

0.53 

49.83 

-1.98 

0.46 

0.98 

70.25 

-3.35 

0.86 

1.82 

AVERAGE  a 

5.27 

8.60 

7.11 

Table  4 —  Smoothed  Pseudorange  Relative  Positioning 


Units  = 
meters 

EAST 

NORTH 

VER1 

1CAL 

MAX 

MIN 

MEAN 

a 

MAX 

MEAN 

0 

MAX 

MEAN 

a 

PIPER 

0.35 

-0.30 

0.08 

0.15 

-0.52 

-0.96 

-0.79 

0.06 

0.99 

-0.53 

0.09 

0.27 

RTE  301 

0.77 

-4.05 

-0.98 

1.26 

3.98 

-1.00 

0.95 

1.11 

3.45 

-6.71 

-1.94 

1.23 

HOLLOMAN 

0.11 

-0.40 

-0.08 

0.16 

0.22 

-0.30 

0.06 

0.11 

0.44 

-0.92 

-0.05 

0.27 

AVERAGE  a 

0.53 

0.43 

0.59 

with  standard  deviations  under  a  meter.  However,  the  smoothed  pseudorange  relative  positioning 
solutions  are  the  best  performers,  with  means  smaller  than  1  m  and  standard  deviations  under  0.5  m. 
This  potential  will  be  explored  in  the  next  section  with  reference  sites  at  greater  distances. 


2.4  Other  Reference  Sites 

PPS  and  PAN  solutions  are  absolute;  that  is,  they  do  not  need  to  have  a  reference  site  collect 
data  simultaneously.  In  effect,  the  reference  sites  for  absolute  positioning  are  those  in  the  GPS 
Monitor  Station  Network.  Since  the  satellite  ephemerides  are  derived  from  observations  obtained  at 
the  reference  sites,  the  satellites  carry  the  information  about  the  monitor  station  locations  and  transfer 
it  to  the  users  through  each  satellite’s  predicted  or  fitted  position.  The  users  then  use  the  satellite 
signals  to  complete  the  link  back  to  the  monitor  stations. 

The  addition  of  simultaneous  observations  at  a  known  static  reference  site  can  be  used  to  re¬ 
move  some  of  the  common  mode  errors  that  exist  in  the  propagation  of  signals  from  orbit  to  the 
ground.  The  uncertainty  of  the  satellite’s  positions  can  also  be  reduced,  thus  forming  a  stronger  link 
to  the  monitor  stations.  During  the  data  collection  by  the  NIMA  vehicle  at  Holloman,  several  other 
static  sites,  some  at  considerable  distances,  simultaneously  collected  data  at  1-s  epochs.  Figure  3 1  is 
a  time  line  showing  the  periods  for  which  data  were  collected  at  seven  other  sites.  The  time  segment 
coinciding  with  the  2-hr  period  for  which  the  F901  provided  kinematic  truth  was  used  to  demonstrate 
high-quality  relative  positioning  over  extended  distances.  The  sites  and  their  approximate  straight-line 
distances  from  Holloman  AFB  are  listed  in  Table  5. 
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The  ellipsoid  heights  of  the 
NIMA  vehicle  were  found  relative  to 
each  of  these  sites  and  are  plotted  in 
the  following  figures.  These  illustrate 
that  relative  positioning  can  give 
acceptable  results  over  continental 
distances.  Figure  32  shows  the  truth 
ellipsoid  height  plotted  with  the  PPS, 

PAN,  and  range-corrected  SPS.  For 
comparison,  the  truth  is  plotted  on 
all  of  the  following  figures.  Figure  33 
shows  the  solutions  with  Austin  or 
Monterey  as  reference  sites.  Fig¬ 
ure  34  shows  the  solutions  with  Cal¬ 
gary  or  Dahlgren  as  reference  sites. 

Figure  35  shows  the  solutions  with 
Ottawa,  Yellowknife,  or  St.  John’s 
as  reference  sites.  The  Canadian  sites 
used  Rogue  receivers,  while  the 
other  sites  used  Ashtech  Z-XIIs.  The  St.  John’s  site  shows  jumps  that  may  indicate  a  receiver 
problem  and  might  have  something  to  do  with  the  10-m  disagreement  for  a  portion  of  the  period 
considered. 

From  the  results  of  this  analysis,  it  can  be  argued  that  the  two-frequency  pseudorange  relative¬ 
positioning  algorithm  is  the  best  all-around  performer.  It  requires  a  reference  site,  but  the  reference 
site  does  not  have  to  be  close  to  the  operational  area;  indeed,  it  can  be  thousands  of  kilometers  distant 
and  still  provide  good  results.  The  benefits  of  relative  positioning  are  that  the  PPS  solutions  are  not 
needed,  and  consequently,  all  receivers  can  be  operated  unkeyed  as  long  as  they  are  capable  of 
providing  two  frequency  pseudorange  observations.  Dual-frequency  observations  are  not  required; 
a  single-frequency  pseudorange  relative  positioning  algorithm  developed  by  Dr.  Remondi  is  used  in 
Section  3. 


f  DATA  SPANS  FROM  ALL  SITES  | 
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Figure  31 — Observation  Spans  from  Other  Static  Reference 

Sites 


Table  5 — Distances  of  Remote  Sites  from  Holloman  AFB,  New  Mexico 


RECEIVER  SITE 

DISTANCE  (km) 

Austin,  Texas  (ARL/UT) 

845 

Monterey,  California  (NPS) 

1476 

Calgary,  Alberta,  Canada  (UC) 

2072 

Dahlgren,  Virginia  (NSWCDD) 

2675 

Ottawa,  Ontario,  Canada  (NRC) 

2926 

Yellowknife,  Northwest  Territory,  Canada 

3282 

St  John’s,  Newfoundland,  Canada 

4622 
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ELLIPSOID  HEIGHTS  OF  NIMA  TRUCK:  ST6 


GPS  TIME  OF  WEEK  (s) 


Figure  32— PPS,  PAN,  and  Range-Corrected  SPS 


ELLIPSOID  HEIGHTS  OF  NIMA  TRUCK:  ST6 
HOLLOMAN  AFB,  8  NOV  96 


Figure  33 — NIMA  Vehicle  Positioned  From  Austin  and  Monterey 
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ELLIPSOID  HEIGHTSOF  NIMA  TRUCK:  ST6 
HOLLOMAN  AFB,  8  NOV  96 


V _ _ _ J 

Figure  34 — NIMA  Vehicle  Positioned  From  Calgary  and  Dahlgren 
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ELLIPSOID  HEIGHTSOF  NIMA  TRUCK:  ST6 
HOLLOMAN  AFB,  8  NOV  96 


V _ _ J 

Figure  35— NIMA  Vehicle  Positioned  From  Ottawa,  Yellowknife, 

and  St.  John’s 
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3.0  NIMA  ROAD  TEST 


During  1 997  on  days  1 90  through  1 93 ,  a  road  test  was  conducted  to  collect  an  extensive  GPS 
data  set  with  multiple  receivers.  The  goal  of  the  4-day  excursion  was  to  collect  GPS  data  while 
traveling  along  interstate  highways  and,  then,  to  evaluate  the  quality  of  the  postprocessed  height 
solutions.  A  NIMA  van  with  three  different  receiver  types  and  two  antennas  traveled  from  St.  Louis, 
Missouri,  to  Toledo,  Ohio,  and  back  along  1-55  and  1-80.  Two  Ashtech  receivers  shared  an  antenna 
mounted  on  the  van  roof  near  the  front  of the  vehicle,  while  a  Trimble  receiver  operated  from  a  single 
antenna  mounted  about  a  meter  behind  the  first.  Data  was  collected  at  1-s  intervals  most  of  the  time, 
except  for  the  last  day,  when  the  interval  was  2  s  on  all  receivers.  At  several  points  along  the  route, 
transects  were  run  on  roads  perpendicular  to  the  interstate  highways.  NIMA  operated  no  reference 
sites  simultaneously  with  the  4-day  trek,  but  some  data  were  collected  at  NSWCDL  and  30-s  data 
has  been  downloaded  from  three  NOAA  CORS  sites  located  near  St.  Louis,  Milwaukee,  and  Detroit. 
These  data  were  used  for  relative  positioning  solutions. 

The  road  test  data  was  characterized  by  many  satellite  dropouts  and  periods  of  poor  signal 
strength  that  result  in  navigation  solutions  that  are  not  as  good  as  they  could  be  under  more  stable 
conditions.  The  data  quality  is  inferior  to  that  from  the  ARL/UT  aircraft  flight  and  NIMA  Holloman 
AFB  truck  experiments  discussed  in  Section  2.  In  order  to  try  and  better  understand  the  reason  for 
the  poor  quality  data  seen  in  the  NSWCDL  Rte.  301  and  NIMA  road  tests,  an  investigation  of 
multipath  was  undertaken.  The  premise  for  this  initiative  was  that  increased  multipath  due  to  vehicle 
motion  might  be  the  cause  of  the  low-quality  navigation  solutions.  It  was  found  however,  contrary 
to  the  authors’  expectations,  that  the  average  multipath  from  the  NIMA  road  test  data  set  did  not 
exceed  the  average  multipath  computed  from  data  collected  at  static  sites.  However,  it  can  be  seen 
from  the  plots  that  follow,  that  the  character  of  the  computed  multipath  indicator  is  qualitatively 
different  between  the  static  and  dynamic  data  sets.  In  particular,  the  multipath  indicator  computed 
from  the  dynamic  data  is  characterized  by  impulses,  while  the  data  from  static  sites  is  not. 


3.1  Performance  of  the  Receivers 

Before  proceeding  with  the  multipath  analysis,  the  performance  of  the  receivers  will  be 
illustrated.  Tables  of  the  number  of  epochs  tracked  on  each  day  and  plots  of  the  number  of  satellites 
tracked  as  a  function  of  GPS  time  follow.  A  single  day,  193,  was  chosen  for  the  plots  because  the 
receiver  performance  on  that  day  was  generally  the  best  of  the  four  days.  The  two  receivers  that 
shared  an  antenna  were  an  Ashtech  MD-XE  single-frequency  receiver,  whose  performance  is  listed 
in  Table  6,  and  a  Z-XII  dual-frequency  receiver,  listed  in  Table  7.  The  Trimble  receiver,  listed  in 
Table  8,  was  a  4000MSGR  dual-frequency  receiver.  As  shown  in  Table  8,  the  Trimble — though  dual¬ 
frequency  capable — was  usually  set  up  to  collect  single-frequency  observations.  The  data  were 
downloaded  when  the  receiver  memory  filled,  so  on  some  days,  multiple  files  were  collected.  These 
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Table  6 — Observation  Epochs  from  the  Ashtech  MD-XII  Receiver 


DAY 

BEGIN 

TOW 

END 

TOW 

NUMBER 

OF 

EPOCHS 

EPOCHS 
WITH  4 
SATS  % 

NUMBER 

OF 

OBSER¬ 

VATIONS 

OBS  OF 

LI 

% 

OBS  OF 
L2 

% 

OBS  OF 
Cl 
% 

OBS  OF 
PI 

% 

OBS  OF 
P2 

% 

190A 

325180 

337627 

12038 

93.6 

63400 

98.7 

98.7 

100 

0 

0 

191 A 

390661 

401890 

10955 

97.1 

65008 

98.8 

98.8 

100 

0 

0 

1 91 B  ! 

402615 

408700 

6022 

97.6 

37390 

99.0 

99.0 

100 

0 

0 

191C 

414526 

424778 

9933 

92.6 

47016 

98.4 

98.4 

100 

0 

0 

191D 

424779 

434138 

8864 

92.3 

50130 

98.5 

98.5 

100 

0 

0 

192A 

489011 

486519 

7319 

91.0 

40283 

98.0 

98.0 

100 

0 

0 

192B 

486818 

494766 

3922 

93.8 

23811 

96.0 

96.0 

100 

0 

0 

192C 

494768 

509176 

6740 

87.6 

35253 

97.3 

97.3 

100 

0 

0 

192D 

509252 

524244 

7117 

87.2 

35739 

96.6 

96.6 

100 

0 

0 

193A 

559712 

582508 

11136 

95.5 

64403 

97.7 

97.7 

100 

0 

0 

Table  7 — Observation  Epochs  from  the  Ashtech  Z-XII  Receiver 


DAY 

BEGIN 

TOW 

END 

TOW 

NUMBER 

OF 

EPOCHS 

EPOCHS 
WITH  4 
SATS  % 

NUMBER 

OF 

OBSER¬ 

VATIONS 

OBS  OF 

LI 

% 

OBS  OF 
L2 

% 

OBS  OF 
Cl 

% 

OBS  OF 
PI 

% 

OBS  OF 

P2 

% 

190A 

325173 

343291 

18089 

98.9 

112732 

97.9 

90.1 

100 

90.2 

90.1 

191A 

397075 

412506 

15433 

99.9 

114780 

98.0 

92.3 

100 

92.8 

92.3 

191B 

414123 

424783 

10651 

99.5 

61788 

96,9 

83.1 

100 

83.6 

83.1 

191C 

424785 

435488 

9445 

100 

60452 

96.9 

85.7 

100 

86.2 

85.7 

192A 

419018 

486661 

7645 

99.8 

53926 

96.6 

78.6 

100 

79.7 

78.6 

192B 

486764 

497367 

5303 

99.9 

39485 

95.9 

83.1 

100 

83.2 

83.1 

||  192C 

499052 

511359 

6155 

99.3 

35561 

94.8 

80.5 

100 

81.0 

80.5 

192D 

511598 

519249 

3827 

99.7 

26504 

94.3 

67.4 

100 

68.6 

67.4 

192E 

519878 

523479 

1802 

99.7 

11627 

97.5 

88.6 

100 

88.9 

88.6 

193A 

559740 

582743 

11489 

99.9 

80790 

95.9 

88.9 

100 

89.3 

88.9 

are  indicated  by  the  letters  A,  B,  C,  etc.,  following  the  day  number  in  the  left-hand  column.  The 
column  headed  “Epochs  with  Four  Satellites”  indicates  the  percentage  of  epochs  where  solutions 
were  possible.  The  columns  headed  by  percentages  of  observations  for  LI,  L2,  Cl,  PI,  and  P2 
jnrliratp  the  number  of  epochs  where  good  data  from  at  least  one  satellite  was  obtained  for  phase  on 
Lj,  for  phase  on  L2,  for  pseudorange  on  L,  C/A  code,  for  pseudorange  on  L j  P  code,  and  for 
pseudorange  on  L2  P  code. 
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Table  8 — Observation  Epochs  from  the  Trimble  Receiver 


DAY 

BEGIN 

TOW 

END 

TOW 

NUMBER 

OF 

EPOCHS 

EPOCHS 
WITH  4 
SATS  % 

NUMBER 

OF 

OBSER¬ 

VATIONS 

OBS  OF 
LI 

% 

OBS  OF 
L2 

% 

OBS  OF 
Cl 

% 

OBS  OF 
PI 

% 

OBS  OF 
P2 

% 

190A 

325062 

343592 

18498 

96.0 

92897 

94.0 

0.1 

99.9 

89.8 

0 

191A 

390884 

424394 

33469 

98.7 

223331 

94.5 

0 

99.9 

90.1 

0 

191B 

427145 

437436 

10288 

97.6 

66330 

90.9 

51.1 

100 

0 

51.1 

192A 

478832 

486376 

7543 

99.3 

52780 

92.0 

43.4 

100 

0 

43.4 

192B 

486830 

511928 

12547 

98.9 

82781 

92.7 

52.7 

100 

0 

52.7 

192C 

511990 

524252 

6085 

96.1 

40443 

92.1 

49.0 

100 

0 

49.0 

193A 

559922 

582638 

11357 

99.4 

79958 

95.1 

57.9 

100 

0 

57.9 

Plots  of  the  satellites  tracked  on  day  193  by  each  receiver  are  shown  in  Figures  36  through 
38.  The  figures  show  the  elevation  angle  for  each  satellite  in  view  during  the  day.  The  vertical  lines 
indicate  where  observations  from  a  satellite  were  not  recorded  or  were  edited.  Figure  36  shows  the 
results  from  the  MD-X13  receiver.  Figure  37  the  results  from  the  Z-XII,  and  Figure  38  the  results 
from  the  Trimble.  The  Z-XII  displays  the  fewest  edited  epochs  while  the  MD-XII  displays  the  most. 
These  results  are  also  typical  for  the  other  days  of  observations. 

Another  indication  of  the  relative  tracking  performance  of  these  receivers  during  the  NIMA 
road  test  is  given  by  the  total  number  of  satellites  tracked  at  each  epoch  and  the  corresponding 
GDOP.  Figures  39  through  41  show  this  information  for  day  193.  The  constant  changes  in  the 
number  of  satellites  tracked  make  it  difficult  for  the  receivers  to  provide  high-quality  observations. 
When  the  number  of  satellites  in  view  is  dependent  upon  the  external  environment,  the  geometric 
strength  suffers,  and  the  solutions  may  be  poorly  determined.  The  GDOP  is  an  indicator  of  the 
solution  quality  due  to  geometry.  Again  it  is  clear,  based  upon  the  number  of  satellites  tracked,  that 
the  Z-XII  performs  the  best  and  the  MD-XII  the  worst. 


3.2  Multipath  Indicator  Model 

The  multipath  indicator  to  be  described  is  a  way  to  quantitatively  evaluate  the  difference 
between  range  difference  observations  obtained  from  pseudoranges  and  the  equivalent  observations 
obtained  from  phases.  These  differences  are  attributed  to  multipath,  though  any  propagation 
phenomena  that  affects  the  two  kinds  of  observations  in  a  manner  similar  to  multipath  may  he  lumped 
into  the  indicator.  The  formulation  that  follows  is  attributed  to  Dr.  Benjamin  Remondi.  In  order  to 
begin,  two  frequency  pseudorange  and  phase  observations  are  required.  Therefore,  use  of  the 
observations  from  the  NIMA  road  test  restricts  the  evaluation  to  the  Ashtech  Z-XII  receiver,  since 
the  other  two  receivers  did  not  record  two  frequency  pseudoranges.  Any  multipath— like  time  delays 
that  affect  the  pseudorange  and  phase  observations  in  the  same  manner  cannot  be  detected  by  this 
indicator.  This  means  that  it  is  possible  to  have  rapid  accelerations  in  the  observations  that  may 
violate  the  tracking  ability  of  the  receiver  tracking  loops,  causing  satellite  mistracking,  without  these 
accelerations  being  detected  by  the  indicator. 
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Figure  36 — Satellites  Tracked  by  the  MD-XII  Indicated  by  Their  Elevation 


Angles 


Figure  37 — Satellites  Tracked  by  the  Z-XII  Indicated  by  Their  Elevation 


Angles 
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Figure  38 — Satellites  Tracked  by  the  Trimble  Indicated  by  Their  Elevation 

Angles 


Figure  39 — Number  of  Satellites  Used  by  the  MD-XII  for  Solutions  and 

GDOP 
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Figure  40 — Number  of  Satellites  Used  by  the  Z-XII  for  Solutions  and 

GDOP 


Figure  41 — Number  of  Satellites  Used  by  the  Trimble  for  Solutions  and 

GDOP 
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The  Li  and  L2  pseudorange  observation  equations  are  given  by  Equations  (1)  and  (2)  below. 


Plji  =  r-ji  +  T-ji  +  -J-ji  +  c(<At"i  ~  At-J +  (1) 

f\ 


P2ji  =  r.Jt  *  T.Jt  +  -  At.jf )  +  M27, 

fi 


(2) 


The  symbols  are  defined  by  the  following: 

yO  is  the  pseudorange  observation  on  L,  to  satellite./  at  time  t, 
r  ..  is  the  geometric  range  to  satellite./  at  time  / 

Tj.  is  the  tropospheric  refraction  to  satellite 7  at  time  t, 

J  is  the  ionospheric  refraction  to  satellite/  at  time  t, 

At  .  is  the  receiver  time  offset  at  time  t, 

••l 

A  t,ji  is  the  satellite  /  time  offset  at  time  tf 

is  the  pseudorange  multipath  at  L,  to  satellite/  at  time  / 


The  difference  between  the  two  observations  is  written  as  Equation  (3).  All  the  geometric  information 
is  removed  in  the  difference  because  both  frequencies  experience  common  delays  due  to  the 
geometric  propagation  time.  The  receiver  and  satellite  clock  offsets  are  also  common  and  are 
removed.  What  remains  are  the  ionospheric  effects  and  multipath. 


Plji  Plji 


\ 


fi 


+  ^2ji  ^ lji 


(3) 


A  similar  set  of  equations  can  be  written  for  the  phase  observations.  These  are  written  as 
Equations  (4)  and  (5),  where  the  integer  cycle  ambiguity  has  been  omitted. 


-  r.„  ♦  T.j,  -  \l.jt  *  -  A,,)  +  ■»„, 

f\ 


(4) 


*1  <P2Ji  =  r-ji  +  T-Ji  ~  4^  +  -  At^  +  ^ 

Ji 


(5) 
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The  symbols  in  these  equations  are  defined  by  the  following: 

is  the  phase  observation  on  L,  to  satellite  j  at  time  t, 
Aj  is  the  Lj  wavelength  in  meters 
r.j.  is  the  geometric  range  to  satellite  j  at  time  /, 

7\.  is  the  tropospheric  refraction  to  satellite  j  at  time  tt 
I.  -  is  the  ionospheric  refraction  to  satellite  j  at  time  t, 

A  is  the  receiver  time  offset  at  time  ti 

A  t'ji  is  the  satellite  j  time  offset  at  time  t, 

m1ji  is  the  phase  multipath  at  Lj  to  satellite  j  at  time  t, 


The  difference  between  the  two  observations  is  written  as  Equation  (6).  The  same  terms  are  removed 
as  in  the  pseudorange  equations  above. 


A24>2ji  Aj  4>lji  -  cl 


■j> 


_j_  _  _1_ 

A' A 


+  m2Ji  -  mXji 


(6) 


Next,  Equations  (3)  and  (6)  are  added  to  remove  the  ionospheric  term.  What  is  left  in 
Equation  (7)  is  just  the  difference  in  multipath  between  frequencies,  and  between  pseudorange  and 
phase.  Usually,  it  is  assumed  that  the  phase  multipath  terms  are  much  smaller  than  those  contaminat¬ 
ing  the  pseudorange,  so  the  indicator  is  taken  to  be  the  multipath  difference  between  the  two- 
frequency  pseudorange.  In  practice,  these  values  at  each  epoch  are  differenced  from  a  reference 
epoch,  so  that  at  the  first  epoch  the  indicated  value  is  zero.  The  deviations  from  zero  are  then 
accumulated  at  later  epochs  as  long  as  phase  lock  is  retained. 

P2ji  ~  Piji  +  A2<t>2ji  ~  A}  <t>iji  =  M2ji  -  MXJi  +  m2ji  -  mXji  (7) 


3.3  Multipath  Indicator  Results 

Two-frequency  pseudorange  and  phase  data  are  required  to  compute  this  multipath  indicator; 
therefore  only  the  results  from  the  Ashtech  Z-XEI  receiver  can  be  reported.  As  described  above,  the 
two-frequency  observations  are  necessary  in  order  to  remove  the  ionospheric  contribution,  and  the 
pseudorange  and  phase  are  used  together  to  eliminate  the  geometric  effects.  What  remains  are  the 
ionosphere-free  differences  between  pseudorange  and  phase.  For  comparison  purposes,  multipath 
profiles  from  several  static  sites  unrelated  to  the  NIMA  data  set  are  presented  first.  Figure  42  shows 
the  multipath  evaluated  at  a  static  site  Patrick  located  near  Frederick,  Maryland,  in  May  1997. 
Several  satellites  are  shown.  Each  curve  starts  at  the  left  at  an  arbitrary  value  of  zero.  Ideally,  the 
pseudorange  and  phase  should  track  each  other,  the  only  differences  being  due  to  noise  in  the 
measurements.  As  time  continues  to  the  right,  differences  are  accumulated  between  the  pseudorange 
and  phase  that  must  be  due  to  the  different  multipath  characteristics  seen  by  the  pseudorange  and 
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phase  observations.  Since  it 
is  generally  agreed  that  the 
phase  multipath  error  is 
much  smaller  than  the  pseu¬ 
dorange  counterpart,  the  var¬ 
iations  seen  in  Figure  42 
must  be  primarily  due  to  the 
difference  between  the  two- 
frequency  pseudorange  mul¬ 
tipath.  Any  multipath  that  is 
common  to  both  data  classes 
and  both  frequencies  cannot 
be  distinguished  by  this  indi¬ 
cator. 

Observations  of  the 
same  satellites  at  the  same 
time  from  a  different  site 
should  show  a  different  sig- 


hplThindSato^plot  fronTa  site  called  Melsage.  Other  sites  were  also  operating  at  the  same  time  as 
Patrick  and  Melsage,  and  all  show  different  multipath  characteristics. 


The  multipath  plot  shown  in  Figure  44  shows  the  signature  from  one  satellite  PRN02,  which 
was  observed  from  several  sites  at  the  same  time.  The  signatures  are  different  at  each  site,  as  wou  d 
be  expected  if  the  variations  are  due  to  multipath.  But  in  all  cases,  the  variations  have  a  quasi-penodic 
nature  around  a  slowly  varying  trend  line.  The  signature  observed  from  the  MIMA  van  contams  spikes 


or  impulses  that  correspond 
to  periods  where  the  satellite 
tracking  was  interrupted. 
The  cause  of  these  impulses 
may  be  due  to  the  receiver 
tracking  loops  trying  to  reac¬ 
quire  weak  or  rapidly  varying 
signals.  The  loss  of  satellites 
is  frequent  in  the  NIMA  van 
data,  especially  on  L2  since  it 
is  inherently  weaker  than  the 
Lj  signal.  Figure  45  illus¬ 
trates  the  impulsive  nature  of 
the  multipath  indicator  for 
day  190,  with  observations 
collected  by  the  Z-XII  re¬ 
ceiver.  On  average,  the  stan¬ 
dard  deviations  of  the  multi- 
path  indicators  for  the  NIMA 
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van  with  impulses  are  about 
the  same  as  for  the  static 
cases.  Table  9  lists  these  sta¬ 
tistics. 


Further  evidence  is 
available  from  the  Route  301 
observations  collected  near 
NSWCDL  in  April  1997, 
which  were  described  in  Sec¬ 
tion  2.  In  this  case,  an  Ash- 
tech  Z-XII  receiver  with  an¬ 
tenna  on  the  roof  of  a  car 
traveled  along  a  four-lane 
divided  highway  collecting 
data  at  1-s  intervals.  The 
multipath  indicators  for  these 
observations  are  shown  in 
Figure  46.  Four  of  the  five 
satellites  were  mostly  contin¬ 
uously  tracked,  while  PRN1 5  was  not.  The  multipath  indicator  for  the  continuously  tracked  satellites 
(PRN  03, 18,  19,  and  31)  resembles  the  results  for  the  static  cases  more  than  the  NIMA  van  plots. 
But  the  result  for  PRN15  is  more  like  the  NIMA  van  results  because  of  the  occasional  impulses, 
which  coincide  with  periods  when  dropouts  occurred. 


Figure  44 — Multipath  Indicator  ofPRN02  from  Five  Sites 


The  conclusion  de¬ 
rived  from  these  data  sets  is 
that  multipath  is  not  the  pri¬ 
mary  cause  of  the  dropouts 
since  the  multipath  indicator 
is  generally  small  for  satel¬ 
lites  that  are  continuously 
tracked,  whether  or  not  the 
site  is  static.  When  satellite 
tracking  is  interrupted  and 
reacquisition  is  attempted, 
the  receiver  may  produce  a 
few  observations  that  cause 
the  multipath  indicator  to 
register  an  impulse.  This  may 
be  especially  true  because  of 
the  lower  signal  level  of  L2. 
The  multipath  indicator  is 
sensitive  to  the  differences 
between  the  two  frequencies, 
or  if  the  receiver  treats  these 


Figure  45 — Multipath  Indicator  for  Day  190  of  the  NIMA  Road 

Test 
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Table  9— Multipath  Indicators:  Static  and  Dynamic  Cases 


|  STATIC  SITES  /  PRN02 

ROAD  TEST  DAY  190 

SITE 

MULTIPATH 

PRN 

MULTIPATH 

PATRICK 

0.11 

06 

0.22 

MELSAGE 

0.19 

10 

0.06 

FAIRFAX 

0.14 

17 

0.06 

JET 

0.13 

23 

0.11 

NBS5 

0.11 

26 

0.17 

signals  differently  because  of  the  differing  signal  strength,  the  multipath  indicator  may  be  indicating 
a  receiver  instrumental  effect  rather  than  a  propagation  effect.  In  either  case,  this  analysis  cannot 
detect  multipath  that  is  common  to  both  data  classes,  so  no  conclusions  can  be  reached  about 
common-mode  multipath  signals. 

When  there  is  solid  tracking,  the  multipath  indicator  produces  small  numbers,  which  indicate 
little  difference  between  the  pseudoranges  and  phase  observations  and  between  frequencies.  An 
example  of  this  occurred  on  a  stretch  of  1-55  near  Lincoln,  Illinois.  Figure  47  shows  the  multipath 
indicator  for  this  segment.  Note  that  because  the  numbers  are  small,  the  scale  has  been  expanded 
relative  to  Figures  4 1  through  46.  These  results  indicate  that  the  key  to  high-quality  solutions  depends 
upon  a  reduction  in  the  num¬ 
ber  of  satellite  dropouts,  so 
that  satellite  geometry  re¬ 
mains  good,  and  reacquisi¬ 
tion  attempts  are  minimized. 

A  practical  implementation 
of  this  may  require  special 
attention  to  detail  in  the  se¬ 
lection  of  the  GPS  antenna 
and  receiver.  Maximization 
of  the  satellite  signal  on  both 
frequencies  implies  that  a 
low-noise,  high-gain  antenna 
preamplifier  with  superior 
out-of-band  signal  rejection 
be  selected  to  eliminate  inter¬ 
ference.  Spatial  diversity 
through  the  use  of  several 
antennas  may  also  be  helpful. 

Figure  46— Multipath  Indicator  for  Route  301  Excursion 
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3.4  Initial  Point 

Before  beginning  the 
NIMA  road  test,  the  van  col¬ 
lected  data  with  the  three 
receivers  while  static  at  the 
North  Annex  in  St.  Louis. 

This  period  of  about  ten  min¬ 
utes,  before  the  van  began  its 
journey,  provides  an  oppor¬ 
tunity  to  view  the  relative 
scatter  in  the  solutions  a- 
mong  the  three  receivers. 

The  Ashtech  Z-XII  naviga¬ 
tion  solutions  are  obtained 
from  two-frequency  iono- 
spherically  corrected  pseu¬ 
doranges.  The  solutions  for 
the  other  two  receivers  are 
single-frequency  solutions 
using  the  broadcast  ionospheric  model.  The  North  Annex  location  was  more  than  40  km  from  the 
CORS  site  STL3 — too  far  to  be  able  to  obtain  On-The-Fly  kinematic  solutions  for  use  as  truth.  Also, 
the  CORS  sites  record  data  at  30-s  intervals  and,  therefore,  has  the  potential  to  provide  only  about 
20  static  truth  solutions  before  movement  begins.  This  is  too  few  to  allow  a  reliable  kinematic 
solution  with  the  current  software. 

The  figures  that  follow  show  the  consistency  among  navigation  solutions  obtained  in  different 
ways.  Figure  48  shows  the  longitude  and  latitude  of  the  simulated  PPS  sohitions  at  North  Annex  for 
each  of  the  receivers.  Figure  49  shows  the  PPS  ellipsoid  height  plotted  against  the  GPS  time  of  week. 
For  reference,  a  1-m  square  is  included  in  the  longitude-latitude  plot.  The  corresponding  PAN 
solutions  are  shown  in  Figures  50  and  51.  Note  that  the  Trimble  solutions  are  separated  from  the 
Ashtech  solutions  by  about  a  meter.  This  corresponds  to  the  separation  of  the  Ashtech  and  Trimble 
antennas.  All  three  PAN  solutions  shown  in  Figure  50  are  shifted  to  the  west  by  about  a  meter 
compared  to  the  PPS  solutions  seen  inFigure48.  The  heights  are  little  changed  as  a  result  of  the  PAN 
improvement.  Since  there  is  no  direct  truth,  the  only  comparison  that  can  be  made  is  through  relative 
positioning  with  respect  the  NOAA  CORS  sites.  The  absolute  position  of  CORS  sites  is  known  to 
be  better  than  about  10  cm  in  all  components. 

Solutions  using  the  data  from  each  receiver,  in  turn  referenced  to  a  CORS  site  and  data 
recorded  at  NSWCDL,  are  shown  below  for  comparison  to  the  PAN  solutions.  Figures  52  and  53 
show  the  MD-XD  solutions.  Figures  54  and  55  the  Z-XII  solutions,  and  Figures  56  and  57  the 
Trimble  solutions.  These  relative  solutions  use  Dr.  Benjamin  Remondi’s  single-frequency  pseudo¬ 
range  relative-positioning  algorithm  with  the  broadcast  ionospheric  model.  The  solutions  are  at  30-s 
intervals  corresponding  to  the  CORS  observation  interval.  Figures  54  and  55  also  mclude  a  range- 
corrected  SPS  solution.  For  this  solution,  the  three  CORS  sites  at  St.  Louis,  Milwaukee,  and  Detroit 
were  used  as  the  weighted  reference  sites  that  generated  the  range  corrections  which  were  applied 


Figure  47— Multipath  Indicator  near  Lincoln,  Illinois 
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NIMA  SRTM ROAD  TEST,  DAY  190 
STARTING  POINT  AT  NORTH  ANNEX  (ST.  LOUIS) 
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Figure  48 — PPS  Solutions  at  North  Annex,  Longitude  vs.  Latitude 
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NIMA  SRTM  ROAD  TEST,  DAY  190 
STARTING  POINT  AT  NORTH  ANNEX  (ST.  LOUIS) 
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Figure  49 — PPS  Solutions  at  North  Annex,  Ellipsoid  Height  vs.  GPS  Time 

of  Week 
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NIMA  SRTM  ROAD  TEST,  DAY  190 
STARTING  POINT  AT  NORTH  ANNEX  (ST.  LOUIS) 
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Figure  50— PAN  Solutions  at  North  Annex,  Longitude  vs.  Latitude 


NIMA  SRTM  ROAD  TEST,  DAY  190 
STARTING  POINT  AT  NORTH  ANNEX  (ST.  LOUIS) 


Figure  51— PAN  Solutions  at  North  Annex,  Ellipsoid  Height  vs.  GPS 

Time  of  Week 
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NIMA  SRTM  ROAD  TEST,  DAY  190 
STARTING  POINT  AT  NORTH  ANNEX  (ST.  LOUIS) 
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Figure  52 — Relative  Longitude  vs.  Latitude  Solutions  Compared  to  the 
PAN  Solutions  for  the  MD-XII  Receiver 
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Figure  53 — Relative  Ellipsoid  Height  vs.  GPS  Time  Solutions  Compared 
to  the  PAN  Solutions  for  the  MD-XII  Receiver 
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NIMA  SRTM  ROAD  TEST,  DAY  190 
STARTING  POINT  AT  NORTH  ANNEX  (ST.  LOUIS) 


V _ J 

Figure  54 — Relative  Longitude  vs.  Latitude  Solutions  Compared  to  the 
PAN  Solutions  for  the  Z-XII  Receiver 
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Figure  55— Relative  Ellipsoid  Height  vs.  GPS  Time  Solutions  Compared 
to  the  PAN  Solutions  for  the  Z-XII  Receiver 
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NIMA  SRTM  ROAD  TEST,  DAY  190 
STARTING  POINT  AT  NORTH  ANNEX  (ST.  LOUIS) 


Figure  56 — Relative  Longitude  vs.  Latitude  Solutions  Compared  to  the 
PAN  Solutions  for  the  Trimble  Receiver 
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Figure  57 — Relative  Ellipsoid  Height  vs.  GPS  Time  Solutions  Compared 
to  the  PAN  Solutions  for  the  Trimble  Receiver 
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to  the  Z-XII  van  observations  at  the  North  Annex.  This  technique  is  much  noisier  than  the  relative 
positioning  results  because  no  smoothing  was  applied  to  the  corrections.  Better  results  can  be 
obtained  with  a  more  sophisticated  algorithm. 

In  all  cases,  the  longitude  and  latitude  relative  solutions  agree  well  with  the  PAN  solutions. 
The  height  solutions  for  all  relative  solutions  are  consistent,  but  the  PAN  height  solutions  indicate  a 
few  meters  lower,  especially  in  the  case  of  the  Trimble  receiver.  Six  satellites  were  tracked  on  the 
three  receivers,  except  for  a  period  near  the  beginning  where  the  MD-XII  receiver  alternated  between 
five  and  six.  This  is  apparent  in  the  GDOP  plot  shown  in  Figure  58.  Since,  in  this  case,  GDOP  cannot 
be  the  explanation  for  the  disagreement  between  the  Trimble  heights  and  the  relative  solutions  or  the 
Ashtech  PAN  solutions,  the  cause  may  be  due  to  using  the  ionospheric  model.  This  is  unlikely, 
though,  since  Trimble  and  MD-XII  solutions  were  single-frequency,  and  the  same  model  was  used 
for  both. 


/ - - 


NIMA  SRTM  ROAD  TEST,  DAY  190 
STARTING  POINT  AT  NORTH  ANNEX  (ST.  LOUIS) 


V _ _ _ J 

Figure  58— GDOP  Computed  While  Static  at  the  North  Annex 


3.5  Transect  1  Solutions 

The  location  for  Transect  1  was  at  the  intersection  of  1-55  with  1-72  near  Springfield,  Illinois, 
as  in  Figure  59.  The  solution  tracks  indicate  the  intersection  to  be  a  cloverleaf,  with  1-55  going  under 
1-72.  There  were  four  passes  made  through  this  interchange.  The  first  pass  was  on  day  190  traveling 
north.  The  second  was  on  day  191  traveling  south.  The  third  was  also  on  day  191  traveling  north.  The 
last  was  on  day  193  traveling  south.  The  passes  on  day  191  were  not  transects— they  were  extra  trips 
to  meet  a  courier  from  St.  Louis  with  a  replacement  download  cable  for  the  Z-XII  receiver.  The 
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Trimble  receiver  recorded  all  of  the 
passes,  the  Z-XII  missed  the  south 
pass  on  day  191,  and  the  MD-XII 
missed  day  190  and  the  south  pass 
on  day  191.  The  solutions  from  the 
MD-XII  receiver  were  of  poor  qual 
ity  and  are  not  included  in  the  fol¬ 
lowing  discussion.  The  solution 
tracks  illustrated  in  the  plots  that 
follow  give  a  qualitative  feel  for  the 
accuracy  of  the  PPS  and  PAN  solu¬ 
tions  obtained  from  the  Z-XII  and 
Trimble  receivers. 

Figure  60  shows  the  Z-XII 
longitude  and  latitude  solutions  at 
this  location,  and  Figure  61  shows 
the  ellipsoid  height  vs.  longitude.  The  same  plots  with  PAN  solutions  from  the  Trimble  receiver  are 
shown  in  Figures  62  and  63.  The  height  solutions  have  enough  resolution  to  show  that  the  east- west 
route  on  1-72  passes  over  the  north-south  route,  though  there  is  some  height  inconsistency  between 
the  solutions  on  days  190  and  193.  Since  a  different  satellite  constellation  was  involved  on  the  two 
days,  the  difference  may  be  GDOP-related.  In  order  to  investigate  this,  GDOP  plots  are  shown  in 
Figures  64  and  65  in  a  format  similar  to  the  height  solutions.  They  indicate  that  there  was  a  lot  of 


Figure  60— Transect  1  Interchange,  Longitude  vs.  Latitude,  Z-XII 

Receiver 
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NIMA  SRTM  TEST,  DAY  190,  191,  193 
PAN  SOLUTIONS,  Z-XII  RECEIVER 
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Figure  61 — Transect  1  Interchange,  Longitude  vs.  Ellipsoid  Height,  Z-XII 

Receiver 


( - ) 

NIMA  SRTM  TEST,  DAYS  1 90,  1 91 ,  1 93 
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Figure  62 — Transect  1  Interchange,  Longitude  vs.  Latitude,  Trimble 

Receiver 
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Figure  63 — Transect  1  Interchange,  Longitude  vs.  Ellipsoid  Height, 

Trimble  Receiver 


Figure  64 — Transect  1  Interchange,  Longitude  vs.  GDOP,  Z-XH  Receiver 
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Figure  65 — Transect  1  Interchange,  Longitude  vs.  GDOP,  Trimble 

Receiver 


satellite  dropout  activity  in  the  vicinity  of  the  interchange  overpass,  with  satellites  being  lost  due  to 
blockage.  The  Trimble  seemed  to  perform  better  by  being  able  to  produce  more  solutions,  but  all 
Trimble  solutions  are  single-frequency.  The  Z-X1I  solutions  were  two-frequency,  and  the  software 
omits  solutions  when  L2  observations  are  lost,  though  L,  may  still  have  been  tracking.  Because  ofthis 
there  may  be  fewer  satellites  used  in  the  Z-XII  solutions;  consequently,  the  Trimble  GDOPs  on  days 
190  and  193  generally  are  the  same  or  better  than  the  Z-XII  solutions. 


3.6  Lincoln,  Illinois 

A  segment  along  155  near 
Lincoln,  Illinois,  is  include  to  show  a 
region  of  good  agreement  between 
PAN  solutions  obtained  on  different 
days.  A  map  of  this  area  is  included 
as  Figure  66.  Figures  67  and  68  are 
PAN  longitude  vs.  latitude  and  longi¬ 
tude  vs.  ellipsoid  height  solutions 
from  the  Z-XII  observations.  This 
part  of  the  interstate  has  some  curves 
with  hills  or  bridges  that  can  be  used 
to  judge  height  repeatability.  Since 


Figure  66 — Lincoln,  Illinois 
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67 — PAN  Longitude  vs.  Latitude  Solutions  Near  Lincoln,  IL,  from 


the  Z-XII 


Figure  68 — PAN  Longitude  vs.  Height  Solutions  Near  Lincoln,  IL,  from 


the  Z-XII 
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tracking  is  good  throughout,  there  must  be  few  obstructions  along  the  route.  There  are  tracks  from 

3  days,  as  in  the  previous  section,  that  can  be  used  for  comparison.  Here,  the  ellipsoid  height 
solutions  show  nice  agreement  among  the  tracks,  with  a  variation  of  only  a  couple  of  meters.  This 
same  segment  was  used  in  Section  3.3  as  an  example  to  illustrate  low-dynamic  multipath. 

Figures  69  and  70  show  the  Trimble  PAN  solutions  and  Figures  71  and  72  show  the  GDOP 
for  this  segment  of  interstate.  These  solutions  are  the  showpieces  of  this  analysis  of  the  NIMA  road 
test.  They  illustrate  that  under  good  conditions  the  ellipsoid  heights  can  be  estimated  to  within  a 
couple  of  meters  with  PAN  improvement. 

The  GDOP  plots  indicate  that,  for  this  segment,  the  Z-XII  has  fewer  cases  where  spikes 
occur,  but  otherwise  the  GDOP  curves  are  similar.  The  day  191  GDOP  curve  located  between  3  and 

4  in  Figure  72  corresponds  to  the  southbound  track  on  that  day;  it  was  not  observed  by  the  Z-XII 
receiver  and  is  not  present  in  Figure  71.  Since  the  results  for  this  segment  are  quite  good,  the 
simulated  PPS  ellipsoid  height  results  are  also  included  in  Figures  73  and  74  for  comparison.  These 
show  a  wider  spread  between  the  days,  at  times  indicating  a  10-m  difference  in  height.  These 
differences  between  PPS  and  PAN  demonstrate  that  the  precise  ephemerides  used  in  the  PAN 
algorithm  do  produce  a  significant  improvement  over  the  PPS  solutions. 


Figure  69 — PAN  Longitude  vs.  Latitude  Solutions  Near  Lincoln,  IL,  from 

the  Trimble 
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re  70 — PAN  Longitude  vs.  Height  Solutions  Near  Lincoln,  IL,  from 

the  Trimble 


Figure  71 — Longitude  vs.  GDOP  Near  Lincoln,  IL,  from  the  Z-XII 
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Figure  73 — PPS  Longitude  vs.  Ellipsoid  Height  Solutions  Near  Lincoln, 


IL,  for  the  Z-XII 
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NIMA  SRTM  TEST,  DAY  190,  191,  193 
PPS  SOLUTIONS,  TRIMBLE  RECEIVER 


V. 


Figure  74 — PPS  Longitude  vs.  Ellipsoid  Height  Solutions  Near  Lincoln, 

IL,  for  the  Trimble 


3.7 1-80  in  Indiana 

The  last  segment  of  the 
NIMA  road  test  to  be  discussed  in 
this  report  is  an  east-west  stretch 
along  1-80  in  Indiana  (see  Figure  75). 
On  day  191  the  NIMA  van  travelled 
east,  and  during  the  return  trip  on 
day  192,  the  van  travelled  west.  For 
the  approximate  25-min  period  cov¬ 
ered,  PPS,  PAN,  and  single-fre¬ 
quency  relative  positioning  solutions 
are  presented.  PAN  solutions  from 
the  Z-XII  receiver  are  illustrated  in 
Figures  76  and  77.  Figure  76  shows 
the  longitude  vs.  latitude  solutions, 
and  Figure  77  shows  the  longitude 
vs.  ellipsoid  height  solutions.  The 


Figure  75 — 1-80  in  Indiana 


vs.  ellipsoid  height  solutions,  i  ne  ,  „„  T,  tra  n  - 

Trimble  solutions  using  L,  plus  the  ionospheric  model  are  shown  in  Figures  78  and  79.  The  tra  g 
If  the  Z-XII  is  better  in  this  case  because  of  fewer  bad  observations which Translate  into  solu tons 
spikes  Also  the  Z-XII  demonstrates  better  consistency  m  the  ellipsoid  height  between  tbeeast  and 
ST  The  offset  between  days  is  on  the  order  of  only  2  to3  m,  whereas  the  Tumble  heigte 
offset  exceeds  5  m  during  a  significant  portion  of  the  segment.  The  GDOP  computed  for  the  two 
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receivers  is  shown  in  Figures  80  and  81.  The  GDOP  for  both  is  at  or  below  four  for  most  of  the 
segment,  exceeding  six  only  during  short  periods.  Notice  that  the  satellite  losses,  in  the  case  of  the 
Z-XII,  may  coincide  with  the  “hills”  in  the  ellipsoid  height  plot.  These  features  must  be  due  to 
overpasses  where  the  interstate  crosses  over  another  road. 

The  PPS  solutions  for  the  Z-XII  are  included  in  Figures  82  and  83,  and  for  the  Trimble  in 
Figure  84  and  85.  These  solutions,  using  the  broadcast  ephemerides,  have  a  10-m  height  offset  in  the 
Z-XII  results  and  are  considerable  noisier,  in  the  case  of  the  Trimble,  than  the  PAN  solutions.  It  is 
questionable  that  height  solutions  with  the  degree  of  variation  shown  here  would  be  useful  for  the 
proposed  application.  PAN  improvement  with  the  precise  ephemerides  seems  to  be  the  preferred 
processing  method  for  the  single  user. 

Single-frequency  relative  positioning  solutions  using  Dr.  Remondis  algorithm  are  presented 
next  These  solutions  are  with  respect  to  CORS  reference  sites  at  St.  Louis,  Milwaukee,  Detroit,  and 
Puerto  Rico,  as  well  as  with  respect  to  the  data  collected  at  NSWCDL.  Figure  86  shows  the  ellipsoid 
heights  on  day  191  from  Z-XII  solutions  compared  to  PAN,  and  Figure  87  shows  the  comparison 
with  the  Trimble  solutions.  Relative  ellipsoid  height  solutions  for  day  192  are  included  as  Figure  88 
for  the  Z-XII  and  Figure  89  for  the  Trimble.  The  solutions  obtained  with  respect  to  the  three  closest 
sites  agree  well  for  both  receivers.  The  results  from  the  most  distant  site— Puerto  Rico— are 
acceptable  on  day  191  but  are  inconsistent  on  day  192.  The  solutions  for  the  Z-XII  receiver  on  day 
192  are  too  low  compared  to  the  rest,  but  they  are  too  high  in  the  case  of  the  solutions  from  the 
Trimble  receiver.  This  discrepancy  is  not  readily  explainable. 
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NIMA  SRTM  TEST,  DAYS  191,  192 
PAN  SOLUTIONS:  Z-XII  RECEIVER 


Figure  77— Z-XII  Longitude  vs.  Ellipsoid  Height  PAN  Solutions  Along 

1-80 


NIMA  SRTM  TEST,  DAYS  191,  192 
PAN  SOLUTIONS:  TRIMBLE  RECEIVER 


Figure  78 — Trimble  Longitude  vs.  Latitude  PAN  Solutions  Along  1-80 


52 


ELLIPSOID  HEIGHT  (m) 


NSWCDD/TR-98/61 


EAST  LONGITUDE  (deg) 


NIMA  SRTM  TEST,  DAYS  191,  192 
GDOP:  Z-XII  RECEIVER 
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NIMA  SRTM  TEST,  DAYS  191,  192 
GDOP:  TRIMBLE  RECEIVER 


Figure  81 — Trimble  GDOP  Along  1-80 


NIMA  SRTM  TEST,  DAYS  191,  192 
PPS  SOLUTIONS:  Z-XII  RECEIVER 


Figure  82 — Z-XII  Longitude  vs.  Latitude  PPS  Solutions  Along  1-80 
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NIMA  SRTM  TEST,  DAYS  191,  192 
PPS  SOLUTIONS:  Z-XII  RECEIVER 


Figure  83 — Z-XII  Longitude  vs  Ellipsoid  Height  PPS  Solutions  Along 

1-80 
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NIMA  SRTM  TEST,  DAYS  191,  192 
PPS  SOLUTIONS:  TRIMBLE  RECEIVER 
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Figure  84 — Trimble  Longitude  vs.  Latitude  PPS  Solutions  Along  1-80 
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NIMA  SRTM  TEST,  DAYS  191,  192 
PPS  SOLUTIONS:  TRIMBLE  RECEIVER 


e  85 — Trimble  Longitude  vs.  Ellipsoid  Height  PPS  Solutions  Along 

1-80 


A 


NIMA  SRTM  TEST,  DAY  191 
RELATIVE  SOLUTIONS:  Z-XII  RECEIVER 


Figure  86— Z-XII  Ellipsoid  Height  Relative  Solutions  on  Day  191 
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NIMA  SRTM  TEST,  DAY  191 
RELATIVE  SOLUTIONS:  TRIMBLE  RECEIVER 


Figure  87— Trimble  Ellipsoid  Height  Relative  Solutions  on  Day  191 
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NIMA  SRTM  TEST,  DAY  192 
RELATIVE  SOLUTIONS:  Z-XII  RECEIVER 
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Figure  88— Z-XII  Ellipsoid  Height  Relative  Solutions  for  Day  192 
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NIMA  SRTM  TEST,  DAY  192 
RELATIVE  SOLUTIONS:  TRIMBLE  RECEIVER 


Figure  89 — Trimble  Ellipsoid  Height  Relative  Solutions  for  Day  192 


4.0  SUMMARY  AND  CONCLUSIONS 

Requirements  exist  for  instruments  that  can  accurately  determine  terrain  heights  when 
operating  from  a  distance.  Proper  evaluation  of  the  accuracy  obtained  from  these  remote  sensing 
devices  requires  a  source  of  ground  truth  over  wide  areas.  In  order  to  support  this  ground  truth 

-  .  ...  a  ,  _ liPlffllt  QOMltlOTlS  IS  OT  intCfCSt.  1  lilS 


The  results  are  very  much  dependent  upon  the  platform  and  the  local  environment.  Satellite 
dropouts  due  to  line-of-sight  blockage  can  be  a  problem  in  the  case  of  ^oimd  velucies^  The  aucraft 
solutions  presented  in  Section  2.1  were  virtually  trouble  free.  For  this  test,  the  PAN  so  utions 
referenced  to  the  kinematic  truth  had  a  mean  vertical  error  of  1.02  m  and  a  standard  deviation 
0  83  m.  The  mean  and  standard  deviations  in  the  other  two  components  were  smaller.  These  should 
be  compared  to  a  mean  of  9.57  m  and  standard  deviation  of  2.90  m  for  the  simulated  PPS. 

The  results  presented  in  Section  2.3  from  low  dynamic  vehicle  motion  over  flat  treeless  terrain 
at  Holloman  AFB  were  also  very  good.  For  this  test  vehicle,  the  PAN 

kinematic  truth  had  a  mean  vertical  error  of -0.76  m  and  a  standard  deviation  of  0.81  m.  The  mean 
and  standard  deviations  in  the  other  two  components  again  were  smaUer.  These  results  should  be 
compared  to  a  mean  of -0.79  m  and  standard  deviation  of  2.01  m  for  the  simulated  PPS. 
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Difficulties  were  experienced  with  data  collected  by  receivers  with  antennas  mounted  on 
vehicle  rooftops  traveling  along  highways.  These  results  were  presented  m  Section  2.2  for  data 
collected  by  NS  WCDD  near  Dahlgren  and  in  Section  3 .0  for  the  extensive  NIMA  cross  county 
test.  In  these  cases  satellite  dropouts  were  common  and  the  quality  of  the  solutions  appeare  o 
dependent  upon  the  receiver  and  antenna  combination.  For  the  Route  301  test,  the  P  so u lions 

referenced  to  the  kinematic  truth  had  a  mean  vertical  error  of -3.19  m  and  a  standard  devia 
5  06  m.  The  mean  and  standard  deviations  in  the  other  two  components  were  on  the  same  order. 
These  should  be  compared  to  a  mean  of -2.80  m  and  standard  deviation  of  5.47  m  for  the  simulated 
PPS.  The  feet  that  there  is  little  difference  between  PAN  and  PPS 

source  of  error  was  not  the  broadcast  ephemerides,  which  can  be  improved  by  the  PAN  technique 
but  some  other  error  source.  The  NIMA  road  test  did  not  have  a  kinematic  truth  &r  compansoi^but 
did  exhibit  many  of  the  same  problems  as  the  Route  301  test.  Day-to-day  repeatabihty  over  the  s^e 
terrain  indicates  that  the  PAN  ellipsoid  height  solution  errors  for  the  road  test  were  at  the  2-5  meter 
level.  The  horizontal  repeatability  was  somewhat  better. 

The  solutions  fromthese  data  sets  do  indicate  that  GPS  can  provide  ellipsoid  height  solutions 
most  of  the  time  while  traveling  along  interstate  highways.  The  accuracy  of  the  height  solntioMcan 
be  quoted  only  in  a  statistical  sense  because  of  the  many  variables  that  influence  the  result.  Adequate 
preparation  of  the  vehicle,  receiver,  and  antenna  are  needed  to  insure  the  best  Perf°™^ 
of  the  system  over  a  known  course  under  various  environmental  conditions,  anff satellite  availa  ty 
should  be  done  to  establish  the  baseline  performance  statistics.  Truth  solutions  from  an  mdependent 
source  are  also  helpful  to  increase  confidence  in  the  baseline  results. 
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